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Abstract 

We analyze the Krawtchouk polynomials K n (x,N,p,q) asymptotically. We use 
singular perturbation methods to analyze them for N — ► oo, with appropriate scalings 
of the two variables x and n. In particular, the WKB method and asymptotic matching 
are used. We obtain asymptotic approximations valid in the whole domain [0, N] x 
[0, iV], involving some special functions. We give numerical examples showing the 
accuracy of our formulas. 

1 Introduction 

Definition 1 Let n, N > be integers. The Krawtchouk polynomials K n {x) are defined by 

m 

where 

< p, q < 1 p + q= 1. (2) 
The binary Krawtchouk polynomials are the special case with p = 1/2 = q J73j /. ITffl . 
* e-mail: dominicd@newpaltz.edu 
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The Krawtchouk polynomials j2TJ are one of the families of classical orthogonal polyno- 
mials of a discrete variable [""""[. They satisfy the orthogonality relation 



k=0 

with weight function 



JV / \ 

Ki(k)Kj(k)g(k)= (,) (pg)^-, i,j = 0,...,N 



n(.r) ----- | ^ )p x q N 



Writing the Krawtchouk polynomials in the extended form K n (x) = K n (x, N,p,q), we 
have the the symmetry formula 

K n (x, N, p, q) = (-l) n K n (N - x, N, q, p). (3) 

They also satisfy the three-term recurrence 

{n + l)K n+1 (x) + pq(N -n+ l)K n ^{x) + [p{N - n) + nq - x] K n {x) = 

which is the main object of our analysis. 

The Krawtchouk polynomials are important in the study of the Hamming scheme of 
classical coding theory [T7], [21], [23, 1221, E3], El- Lloyd's theorem [23] states that if 
a perfect code exists in the Hamming metric, then the Krawtchouk polynomial must have 
integral zeros [Sj, [5], [23! • Not surprisingly, these zeros have been the subject of extensive 
research gj, [7j 0, [JU], [II], 02], [US]. 

The Krawtchouk polynomials also have applications in probability theory [S], queueing 
models j^j, stochastic processes EI], quantum mechanics [2], |21] EH] and biology [T^j . 

The asymptotic behavior of the Krawtchouk polynomials as iV — > oo was studied by 
Sharapudinov for x ~ Np and n = O (iV 1 / 3 ) in E2J. He derived an approximation in 
terms of the Hermite polynomials (see also EH] for a similar formula). The general case 
with x,n = O(N) was investigated by Ismail and Simeonov [T2 3 and a uniform asymptotic 
expansion was derived by Li and Wong [22], both using the saddle point method. 

The purpose of this paper is to take a different approach, based on the recurrence formula 
that the Krawtchouk polynomials satisfy and using singular perturbation techniques [26J, 
[3*7] to analyze it. We scale x = yN, n = zN and obtain asymptotic approximations to 
K n (x) for (y, z) G [0, 1] x [0, 1]. Our results agree and extend those obtained in [12], [22! an d 

EH- 

In Section 2 we review the basic properties of the Krawtchouk polynomials. In Sections 3- 
10 we obtain asymptotic expansions from the recurrence formula by using the WKB method. 
We must consider twelve relevant regions of the two-dimensional state space. In Section 11 we 
summarize our results and numerically compare our approximations with the exact formula. 
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2 The WKB approximation 

To analyze the recurrence 

(n + l)K n+ x(x) + pq(N -n+ l)K n -i{x) + \p(N -n)+nq-x] K n (x) = (4) 
subject to the boundary conditions 

K {x) = 1 (5) 

W*)=(/ +1 ) (6) 
K n (0) = (^j {-pf (7) 
K n {N)=( N \ n (8) 



n t 

for large AT, we introduce the scaled variables y, z defined by 

x = yN, n = zN, < y, z < 1. (9) 
We define the function F(y, z) and the small parameter e by 

e = -i AT n (x) = F(ex,en) = F(y,z) (10) 

and observe that i<T n ±i(a;) = z ± e). 
Substituting (|gjl- lfTTTjl in (@J) we get 

(z + e)F(j/, z + e) + pg(l - z + e)F(y, z - e) + [p(l - z) + zg - y] F(y, z) = 0. (11) 

To find F(y, z) for e small, we shall use the WKB method [2E]- Thus, we consider solutions 
which have the asymptotic form 

F(y, z) ~ e v exp [e'^y, z)} L(y, z). (12) 

Using (P| in (fTITl . with 



e V(y, 2 ± e) = £ V(l/» ± $z(v, z) + 2^ M (l/, z)e + O (e 2 ) , 
dividing by exp [e~ l ip(y, z)\ and expanding in powers of e we obtain the eikonal equation 

zU 2 + [p-y + z(q-p)]U +pq(l- z) = (13) 
and the transport equation 

[zU 2 - pq(l -z)]L z + S^ [zU 2 + pq(l - z)] ip zz + U 2 +pqj L — (14) 

with 

U(y, z) = exp [ip z (y,z)] . (15) 



2.1 The functions ip and L 
From (IT5l) we have 



1>(y,z) = Jhx [U(y,z)]dz = z\n(U)- J ' z^dz = z\n{U) - j ^j^jdU 



and using (JT3j) we get 



Uz dz = f (y-p)U-pq 
Z UdU J (U + q)(U-p) 



dU 



1 , y - 1 y 



U U U + q 



dU. 



Hence, 



^(y, z) = In [C/^ 1 (U - p)^ y (U + qf] + A{y) 



where the function A(y) is still unknown. 
From (|14j) we have 



L(y, z) = B(y) exp 



1 [zU 2 + pq(l - z)\ ijj z2 + U 2 +pq 



dz 



zU 2 — pq(l — z) 

and from (|15j) ip zz = U z /U. After changing variables from z to U, we obtain 



L{y,z) = B(y)exp 



(p-y)U + pq 



_2U U + q (p - y) U 2 + 2pqU + pq(q - y) _ 



dU 



Thus, 



L(y,z) = B(y) (U + q) 



U 



(p — y)U 2 + 2pqU + pq(q — y) 



where B(y) is to be determined. 



2.2 The function U 

Rewriting as 



U 2 +(^ + q-p)U+(U f = 



and solving for U we get 



where 



Un(z) 



pq(l - z) 



The discriminant in (|18|) vanishes if ^—^ + q — p = ±2C/q, which is equivalent to y = Y , 
with 

Y ± (z) = p + (q - p) z ± 2zU . (20) 
Rewriting the equation ( 2-^ + g — pj — 4 (Z7q) = as 



1\ 2 / l\ 2 , , f 1\ / 1 



2J + r~2j + 2 (p-^) ^- 2J v z " 2J " w = (21) 

we can see that 

{(y,z) G [0,1] x [0,1] : 2/ = F ± (^)} = E (22) 

where E is an ellipse centered at (|, |) (see Figure |IJ). After rotation by ±? and translation 
to the origin, E reduces to one of the canonical forms 

y 2 z 2 

p<q 

-2 -2 1 1 

y +z = -, p = ~ = q 

y 2 z 2 

H 77T = 1) P > Q- 



p/2 q/2 

The ellipse E is contained in the square [0, 1] x [0, 1] and intersects the y and z axis at the 
points (0,p), (1, q), (p, 0) and (q, 1). Its left side between the points (q, 1) and (p, 0) coincides 
with the curve Y~(z) and its right side with the curve Y + {z). 

For points (y, z) located outside E, U ± (y, z) are real and for points (y, z) inside E U ± (y, z) 
are complex conjugates. When y = Y ± the two values U coalesce and we have 

U + (Y + ,z) = U = U-(Y+,z), U+(Y~,z) = -U = U~(Y-,z). (23) 

Writing the function L(y, z) in terms of Uq we have from (fTTj) 



L(y, Z )=B( y ) X U ~ T ^\;^ - (24) 



\ U-p)(U + q) 
z {U 2 - VI 

In the rest of this paper we shall use the following notation 

^(y,z) =\n\(U ± Y~ 1 - pY' y + q) y ] (25) 



and 



L {y ' z) -i z[(u±f-us] • (26) 

Hence, we write 

K n (x) ~ e"B-(y)exp [s^iy, z) + e~ l A~(y)] L~(y,z) (27) 

+ e"B + (y) exp [e'^+iy, z) + e^A+iy)] L + (y, z) 

where B^iy) and A^iy) are functions to be determined. From (J26|) we see that L (y, z) are 
singular when = Uq, i.e., for y = Y^, and also for z — 0. Therefore, we need to find 
asymptotic solutions valid in those regions. 
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Figure 1: A sketch of the ellipse E and the curves Y (z). 
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3 The boundary n = (Region I) 

For n = 0(1) we have from (JIJ that K n (x) = O (N n ) as N — > oo. Thus, we introduce the 
function Rn (y) and consider solutions of which have the asymptotic form 

K n (x) = N n R^ (*) . (28) 
Using (}2*H|) in (£Q) and expanding in powers of N gives, to leading order 

(n + l)R% + (p-y)RW = 0. 
Solving the recursion above with the initial condition (jSJ) we get 



ni 

and hence 



K n (x) ~ J#)(y) = N ni ^-P)l, n = 0(1). 

n! 



(29) 



3.1 The corner layer at (p,0) (Region II) 

We shall now find an asymptotic solution in the neighborhood of the point (p, 0). We intro- 

(2) 

duce the stretched variable rj and the function Rn (77) defined by 



V =p + r]V 2 P ( l £ V = 0(1), (30) 

Using (J5U|) in (HJ) yields, to leading order, the equation 

R% - 2 V RW + 2n J Ri 2 2 1 = 
which we recognize as the recurrence relation for the Hermite polynomials. Thus, 

Kn(x) ~ Ji^fo) = ^ (f)" 72 ^ (77) (31) 

for n = 0(1) and y — p = O (z 1 ^ 2 ) , where if n (77) is the Hermite polynomial of degree n. 

4 The lower corners (Regions III and IV) 

Setting n = z/e in ()29|) and letting e^Owe obtain 

Jtf> (y) ~ £ 1/2 ^= exp {s- 1 [1 - ln(z) + ln(y - p)] z} (32) 



where we have used Stirling's formula [1 



T(x) ~ J—x x e- x , x^oo. (33) 



From (J2H)-(j2ni) we have as z — > 



[1 — ln(2;) + ln(y — p)]z, y < p 

(l-^ln^+.ln^+lnQ^)*, »>„ < 34 > 

L = 4 . ^/ \ ( 35 ) 



y-p 



+ 0(z), y>p 



and 



[1 — ln(z) + ln(y — p)]z, y > p 

L^J^ + OU), y< P , (37) 
[ z~ 1/2 + 0(z), y>p 

Matching (^ - ([37)1 and flH2} we conclude that 

\K + (^), Y+{z)<y<\ ^> 



with 



K-(y, z) = £ 1/2 -L exp [e^iy, *)] z) (39) 

V Z7T 



*T + (y,z) = e 1/2 ^=exp [e-V + (j/^)] £ + (^)- (40) 



Remark 2 In the remainder of the paper, we will find asymptotic formulas only in the region 
< y < Y + (z), < z < 1. T/ie corresponding results for Y + (z) < y < 1 can &e obtained 
by using the symmetry formula J3J) and noting that under the transformations y — > 1 — y, 
p ^ q, K n — > (— l) n K n , we obtain 

— > ^ + + Z7ri, L~ — > L + 
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5 The boundary x = 



We shall now consider the case x = 0(1). 
Lemma 3 Let x = m, with m an integer, m ^ N. 
1. 



K n (m) = (-p) n (^j (l - ^N- 1 ^ m , m = 
2. Ifn = O(N), then 



Proof. From (JJJ we have for x = m integer 

/ m\ (N — m\ , k 



and (|41|1 follows for m = 0, 1. 

Setting n = ziV and using we get 



/ N—m\ , s k 

Using the above in pHj) we have 

/ T> T\ 771 / \ /• 



«-)~(-prOo-.ri:(T)(-T 



-p)" ( N ) (i - *>* 



n 



k=0 

-i m 

p — Z 



z p 



p(l-z) 



and (j42j) follows. ■ 

Using (J33|) we have, as iV — ► oo 



: " P) '' l - ) ~ V2^(l-^ eXP[£ " V ° W1 



where 

0o(z) = [z — 1) ln(l — z) — z Yn.(z) + z ln(— p). 
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From (jHU) we get, as y — > 



K-(y,z) 



,1/2 



exp 



J0 (z) + \n[^)y 



, < z < 



V^\/z(l - z) 

Hence, as y — > 0, K~(y, z) satisfies the boundary condition (jJJ) for < z < p. 

5.1 The boundary layer at x = 0, p < z < 1 (Region V) 
Taking (jj2j) into account, we define the function i?„ (x) by 



KJx) 



N 



n 



Using (g7D in (gj) yields 

p{n - N)R% + {pN + (q-p)n- x] - nqR { ^_ x = 0. 
Writing in terms of z, e and the function (x, z) defined by 

Rf>(x) = G {5 \x,en) 

we have 

p{\ - z)G (5) (x, z + e) + [z (p - q) + xe - p] G {5) (x, z) + zqG {5) (x, z - e) 
Using the WKB anszat 



G (5 \ 



x, z) ~ e exp 



~<p\x,z) 



W{x,z) 



in (jSUJ), we obtain the equations 

[exp {(j) z ) - 1] [p(z - 1) exp (0 Z ) + zq] W = 

and 

[zq + p{z - 1) exp (20,)] W z 

+ { \ ^ Z ~ ^ 6XP ~ ^ ^ zz ~ X eXp ^ } W = °" 
Solving (j52 > ([53 |) we get 

G (5) (x,z) ~e Tl A (5) (a;)exp [e _ Vi(aO] 
+e T2 J B (5) (x)exp [e - V2O, z)] (z - p)^ 1 y/z(l-z) 
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where 

4>2{ x -, z ) — (1 ~ z ) In (1 — z) + 2 \n(z) + [z — 1) ln(g) — z ln(p) + <p2(x) (55) 

and the coefficients t\,t% and the functions A^\x), B^\x), <fi(x), tp 2 (x) are to be deter- 
mined. 

Using (jUj) we obtain 

P l/2+n 



K n {x)~K® (x,z) 



2tx^/z(\ - z) 



exp 



<j>o(z) + (pi(x) 



A(x) (z — pf 



(56) 



s 



1/2+T2 



exp 



(z - 1) ln(g) + itiz + (p 2 (x) 



B(x) (z — p) 



-sc-1 



5.2 The corner layer at (0,p) (Region VI) 
For x ~ and n w 7Vp, we scale n as 

n = iVp — u\/pqN, u = 0(1) 

and introduce the function G^(x, u) defined by 



KJx) 



N ) {-pTg^ (x Np ~ n ^ 



VpqN ) 



(57) 



(5? 



Using (JH7|) - (JH%|) in (|3%J) we get, to leading order 

d 2 G (e) qq(6) 



du 1 



du 



+ X G^ = 0. 



Solving (|59j) we obtain 



G^ 6 \x, u) = exp 



-u 



(59) 



(60) 



where D x {u) is the parabolic cylinder function and A^(x), B™(i) are functions to be 
determined. Since 

Np-n p-z 



u 



VpqN v / pq£ 

we note that the limit u — > oo corresponds to the matching between regions VI and III, 
while the limit u — > — oo corresponds to the matching between regions VI and V. As z — > p 
we have 



-pY [ N n >~ 



exp 



u 
"~2 



x exp 



7rip - gin (g) i 7 ™ + m 



(62) 
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Thus, 



K n {x) ~ K {6 \x,u) 



V^npq 



[A^\x)B x (u) + B^(x)D x ( 



x exp 

Using US]} in yields 



nip — gln(g) u^/pq [id + In (q)] u 1 
e y/e 4 



K- (y,z) 



\/2npq 




u 



x exp 



7rrp-gln(g) Uy/pq[n\ + \n.(q)} 



u 



Using the well known asymptotic approximation 



u 



Dji(m) ~ exp ( — — ) u x , u — > oo 



in (J53)l and comparing with we conclude that 

A {& \x 

and B^\x) = 0. Therefore, 

K {6 \x,u 

x exp 




,1/2 




DM 



^/2npq 

nip — q In (g) u^/pq [ni + In (g)] 



u 



for x = 0(l),z-p = (e 1/2 ) . 
Using the formula 



D x (— u) ~ exp 
in (jMj) we have 



it 



u cos 7T3; 



2 _ _ /u 2 

xr(x) sin(7rx)w x 1 exp ( — 



x exp 



.1/2 



V2npq 



u 




COS 7TX 



7rip — g In (g) ^v/pg [?ri + In (g)] w 2 

Y 



1 e 1 ^ 



1 /gF 
uV p 



xY(x) sin(7rx) 



x exp 



u ny/pq 

nip — g In (g) u^/pq [ni + In (g)] 
e s/e 



u — > oo 
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Using (jHTjl in (fo7)j) gives, as z { p 



x exp 



Kip — q In (g) + (x) w ^/pg [vri + In (g)] 



it 



+ exp 



£ 2 
7rip — gin (g) + <^ 2 (x) w v /pg [7ri + In (g)]' 



x ■ 



£ 

1 £ T2 



:S (5) (x) (uy/pqe)' 



u ^J2npq 

Matching (JHEJ) and dHHJ) yields n = 0, tp x (x) = 0, r 2 = |, <p 2 (a:) = and 

A (5) (x) = p _a: cos(7rx) 



Thus, 



K (5) (x, . 



,1/2 



2tt v /,z(1 - z) 



COS 7TX 



r- 










"0o(z)~ 


1 exp 









£ X 



K Z — p 



-T(x) sin(7rx) 



gs 



z — p 



exp 



2 — p 
p 

(z — 1) ln(g) + kiz 



for x = 0(1) and p < z < 1. 

From (fTOj) we have, as x — > oo 

if (5) (x, z) - 



,1/2 



2 ^ 
—e\ sm(7rx 

V 7T — p 

In terms of y = xe, (J7TJ) reads 



271-^(1 - z) 
qex 



COS! KX I 



2! — p 
P 



exp 



0o(*) 



z — p 



exp 



— 1) ln(g) + kiz 



x 



if (5) (x, z) 



,1/2 



Ky 

cos I — I exp 



27^(1 - z) V e 



E / Z — p 



exp 



> - 1) ln(g) + vri2; + yln f^LJ - y 
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6 The left upper corner (Region VII) 

We now consider the region <C y < Y~(z), p < z < 1. From (j25j) - ()2f)j) we have for 
y -> 0, z > p 



*P + (y, z) ~ O («) +2/ ln 



2 — p 



p 



my, 



■0 (y, 2) ~ (jz — 1) ln(g) + niz + y In 



qy 

z — p 



y - my, 



and 



L + (y,z) 



L (y, 2;) — i. 



Using (J?3*j) - (J7Hj) in (|27|) and matching with (JZ2J), we conclude that 



K n (x) ~ if (7) (^) = exp ( ^ j [cos (^J K+(y,z) + 2isin (^J K~{y,z) 



T\y 



ixy 



(73) 
(74) 

(75) 
(76) 



which can be written as 



K {7 \y,z) = l - 



cxp 



271% 



K + {y,z) + 



cxp 



27riy 



K-{y,z). (77) 



7 The transition layer y = Y (z) 



As we noted before, the functions ^{y, z) are infinite on the curves y = Y ± (z). Hence, we 
need to find transition layer solutions there. 

7.1 The lower part < z < p (Region VIII) 

We introduce the stretched variable j3, defined by 

y = y~(z) -(3e 2/3 (3 = 0(1), 0<z<p. (78) 
Using (fTSj) in (|23|) - (|2^|) and expanding in powers of e, with (3 > 0, we get 

2 fU Q 1 



i) (y,z) ~ i/j (z) + In 



Uo+P j ^ 
U -q 




3 V z (U +p)(U -q) 



(79) 



where 



and 



^o(z) = «7ri + (z - 1) In (C/ ) + Y~{z) In (C/ - g) + [l - F~(z)] In (U + p) (80) 

L~(y,z) ~ -L^-V4^ (t/o +p) (f/o _ g) ([/o) -3/4 / j-l/4 e - 1 /6 i (gl) 
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Hence, 



K-(y, z) ~ e^V^o + P) Wo ~ Q) (^o)~ 3/4 /T 1/4 



(82) 



x exp 



e Vo(X) + m 



-1/3 



2 /t/i 








3 V z (U +p)(U -q) 



03/2 



for y ] Y (z). Thus, we introduce the function G^ 8 \f3, z) and consider solutions of the form 



K n (x) = e us exp 



ipa(z)£ 1 + ln 



Uq+p 



(y~ - ex) e- 1 



Using (|83j) in (jlj) and expanding in powers of e we obtain 

,(8) _ 2/3 



G® 



Y — ex 



r2/3 



, en 



13) 



G 



where 



Solving we get 



[(y-)T[(p- ? )^ + y--p] 

Q(z) = 



g^ = pe 2 c^ 



Un 




z (U +p)(U -q) 



G®(J3,z) = A®(z)M [Q 2 / 3 f3] + B®(z)Bi [Q 2 ' 3 f5] 



»4) 



*5) 



56) 



where v4( 8 )(2:), 5^ 8 ^(z) are functions to be determined and Ai (•) , Bi (•) are the Airy functions. 
Using the formulas pQ 



Ai (x) 
Bi (x) 



2^ 
1 



/: 1/4 exp ( -^ x3/2 



x — > oo 



(87) 



-1/4 



7T 



x ' exp l -x 



2 3/2 



x — > oo 



in (}86|) and matching with (}82j) we conclude that u 8 = 1/3, i?®(z) = and 



(Uo+P) (Uo-q) 
z(U f 



V3 q-1/3 



Therefore, 
for y-Y~(z) = 0{e 2 / 3 ) 



z) = e 1 / 3 exp 



e 1 ip Q {z) + ln 



Uo+P 
Uo-q 



f3e 



-1/3 



Ai [Q 2 / 3 (3] 



0-1/3 



(89) 
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7.2 The upper part p < z < 1 (Region IX) 

Using (f78j) in (}25j) we have, as y | ^"(2) 



'0 (y, 2;) ~ ^0(2) + In 



Uo+P 
U -q 



*p + (y,z) ~ Vo(^) + 
and from ()26|) we have 



\U -qJ 3 



i - ( -^717ib ir,/VI/6 - 

L+(y, 2 )~4^= ; r 1/ V 1 " i 



where 

#(z) = _e(s) 

Using (jnOD-dHSD in JZZJ we have 



^0 




1 



z (u + P )(q-u y 



cxp 



as ?/ | V (-2), where 



2v^F v^CTqW 



4>o(z)e 1 + In 



U -q 



Pe 



-1/3 



-A + (/?, z) exp ( — -i!}[3 3 / 2 J + iA z) exp 



A ± (/3, z) = exp 



2vri [Y- (z) - /?e 2 / 3 ] 



± 1. 



We consider the ansatz 

K n (x) ~ K {9) ((3,z) = e ua exp 



e Vo(^) + In 



Uo+P 
Uo-q 



Pe 



-1/3 



x [\+(p,z)A^(z)Ai {d 2/3 p) + A-(/3,2)5( 9 )(^)Bi (tf 2 / 3 /?)] 

for y rs p < z < 1 and unknown functions A^(z), B^\z). Using (|87j) 

have, as (3 — > 00 



^ (9) (A 



~ e J exp 



e Vo(z) + In 



Uo-q 



0e 



-1/3 



x 



A + (/3,z)A( 9 )(z)-V^- 1 / 6 r 1/4 exp (-^P 3/2 ) 



A" (/?, z)5 (9) (z) --Ltf-^p- 1 / 4 exp 



+ 



'-$P 3 ' 2 
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Matching (jHHJ) and (JHSJ) we get 



^9 



1 

3' 



^ 9 >(z) 



1 0-V3 

2~7^o' 



i?( 9 )(; 



1 i?-V3 
-i- 



2 



o 



Thus, 



cxp 
2~7^o 



.-Vo(.)+ln(^)/?.- ] S 



x 



[A + (/3, z)Ai ($ 2/3 f3) + VC(J3, z)Bi ($ 2/3 [3)] 



(99) 



8 The interior of E (Region X) 

We shall now find an asymptotic solution for Y~(z) < y < Y + (z). We set (3 = —(3, (3 > in 
(JZSIl-dHD and obtain, for < z < p 



and 



i> (y, z) ~ ?A 



L {y,z) 



*{z) - In (^±P) frfi + W 2 . 
\Uo-qJ 3 



Similarly, from (|25j) - (j2fij) we obtain, for < z < p 



e 4 



1^-1/4^-1/6^ 



\U -qJ 3 



and 



L + (y,z) 



TTi/O-l/4 -1/6 



Introducing the function K^ 10 \y, z) defined by 

K^°\y,z) = e^A^\y)exp [e^-faz)] L~(y,z) 
+e^B ( - w \y)exp [e'^+iy, z)] L + (y,z) 

with A^ 10 \y), B^ 10 \y) to be determined, we have 



(100) 



(101) 



(102) 



(103) 



(104) 



#(io) ( 



y,z) ~ e 



-10-1/6 i 1 £-1/4 exp 



t/o-g 



0e 



-1/3 



(105) 



^4 (10) (l/) exp ( 2 -&^ 2 - ~7ri) + B^(y) exp ("jj©i/3 3/2 + ^ri 
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From (J89J) we have 

AT (9) CM ~e 



^- l_^-V4 ^ 



x sin Qe^ 2 + ^ 
where we have used the asymptotic formula jTj 



e Vo 0) - In 



^ ) 'Be- 



(106) 



Ai f-x) ~ 



-1/4 



7T 



sin I -x 3/2 + 



7T 



X — > OO. 



(107) 



Matching ()105|) and (|106|) we get z/ 10 = \ and 



1 

71 



A^\y) exp ( ^6i/5 3 / 2 - + B^ 10 \y) exp f-^6i/3 3 / 2 



^sinf-6/3 3 / 2 + -, 



exp 



from which we conclude that 



A^\y) 



|0i/3 3 / 2 - hi) + exp f-|ei/? 3 / 2 + ivri 



5 (10) (Z/) 



2,7T 



and therefore 



- ^ 10 )(y, z) = K+(y, z) + z) 



(108) 



for F-(^) < y < Y + (z). 



8.1 Matching the interior of E and the upper part of the transition 
layer Y 

We shall now verify the matching between ()108|) and (}9*9"j) . For p < z < 1 we have 



< ■-(//. : i - (1- i /fe 2 / 3 - *) 27 ri + V-b(s) - In ( (3e^ - ~tfi/3 3 / 2 e, 

U — qj 6 



L (y,z) 



1 1 



V2^/zU^d 
^ + {y,z) ~ ^o(^) - In 



e ^-i/4 e -i/e 



Uo-qJ 3 



and 



1 1 
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-1^-1/4^-1/6^ 



(109) 

(110) 
(111) 

(112) 



Therefore, 



U — q 
2 



\ + (-(3,z) + \-(-/3,z) ( 2 ~ 3/2 1 \ /2 Q ~ 3/2 1 

— exp --$i/3 3/2 + -7ri + exp -$i/3 3/2 - -tti 



Using (HOZD and P 



-1/4 



Bi ( — x) ~ — =- cos I -x 3 / 2 H — | . ./• -• oc 
/tt V 3 4 



in (|9*9*j) we have 



exp 



£ Vo(^) - In 



Uo+P \ ;,.__i /3 



X 



1 & 



-1/4 



2v^F v^ZTq^ 



A + ( A :)sn, ( ^/3 3 / 2 + | j +iA-(-/3,z)cos + | 



or 



z) ~ e ' exp 



£ ^o(^) - In 



1 (3-^ 



A + + A~ ( 2 a -o3/2 1 A A+-A- /2 ~ 3/2 1 . 

exp --$i/? 3/i + -vri H exp ( -$i/T /2 - 7 7ri 



l 

4 



A+- A" 



(113) 



(114) 



2 V 3 

Since from (|96|) we have 

we see that ()114|) agrees with (jllHj) . 

8.2 Matching the interior of E and the corner layer at (p, 0) 

Using (|33p and the asymptotic formula jTHj 

H n {x) ~ v^V^exp - |) cos (v^x - ^ 
in (J3T|) we have, as n — > oo 



n7r\ 

- J , n — > oo 



I r 77 n 

exp < — + - 

n7T 2 2 



£n 



. / — nir 
cos ( v 2nr/ — 



(115) 
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Using (|3Uj) in (|23j). with z = ne we get 



1 ± n 
£ 2 L 



1 + In ( — ) ± 7ri 

en 



^ v2w/i + 



Similarly from (j2b|) we find 



2n 



Using (fTTfi |) -(fTT7 j) in (fTUHj) we obtain 



1 ( T] n 



which is in agreement with (jf f 5j) . 



1 + ln 



pq 



en 



} cos (^/2nrj — j 



9 The boundary layer at z = 1 (Region XI) 

We now consider solutions of (@J) with n ^ N. We introduce the variable j and 
Tj\y) defined by 

n = N-j, jeZ, J>-1, 



Using (jll9|) in (jlj) we have, to leading order, 



qTfH -(q-y) Tf > = 



which can be solved to obtain 



where (y) is a function to be determined. Thus, we have 

K n (x) ~ A(") (y) (-p)^ (l - ^ 

We now introduce the function Rj(y) defined by 



and note that from (??) we have 



i£?(v) = I- 
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Using ()118|) in (jU we obtain, to leading order, 



(q-y)R^l-(l-y)R 



(11) 



which together with ()123j) implies 



i-y 



i+i 



y 7^ ?• 



Therefore, 



Ny\ fl-y 



n J \q — y 



3+1 



yftq. 



From ()122j) and ()125j) we get, for y ft q, 



K n {x) ~ ^(y) = AW (y) ( ^ (-p)^ f 1 - *Y 4 



Ny \ fl-y 



N - j J \q-y 



3+1 



(124) 



(125) 



(126) 



9.1 Matching with the upper corners 

We shall now determine the function (y) by matching (|126|) with (J76j) . As N — ► oo we 



get 



Aj u) (y) ~ AM (y) -±=j exp j^ 1 In (-p) + 



I - In (./:) - luf -//) In ( 1 - ^ 



(127) 



f y ln(y) + (1 - y) In (1 - y) + 7ri 



. i/2 [z V( 1 -y)y . (AT \ 

+e 1 \ — — walNiry) exp 

\ TV y-q 

Using (|118|) in (J25J) we have as £ — > 0, < y < q 

1 - In (je) - ln(-p) + In ( 1 



J In (y - q) 



ip + (y, z) ~ ln(-p) + y In y- j - + 



je (128) 



i> (y, z ) ~y My) + ( 1 ~y) ln (i - y) + (i - y)^ - i ln (y - ( 129 ) 

and from (J26|) we get 

£ + (y, z) ~ -^U z ) ~ s/y{i-y)—i. (130) 

vj^ 9 - y 

Using ipBjl -lfrcn ) in (fZEJ) we have 



y,z) ~ cos l — 



7ry\ 1 



ln(— p) + y ln 



exp 



1 - ln (je) - ln(-p) + ln [ 1 - - 

q 



(131) 



+v / y(i-y) 



i 

y-9 



sin 



ny\ 2e 



e / V 7r 



exp 



y ln(y) + (1 - y) In (1 - y) + vri 



j ln(y - g) 
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Matching (|127|) and ll31l we obtain 



A^\y) 



cxp 



COS 



Hence, we conclude that 



N ) (-pf-> ( q - 



for z = 1 - 0(e) and y ^ q. 



9-2/ 



(132) 



(133) 



9.2 The corner layer at (g, 1) (Region XII) 

The approximation ()133|) ceases to be valid for y m q. Therefore, we need to find another 
expression, which holds in a neighborhood of the point (q, 1). We introduce the new variable 
£ defined by 

y = q + £^ql, £ = 0(1). (134) 
Using ()134|) in ()133|) we have, as e — > 



(y) ~ ex P 



pln(p) + gln(g) /2po /g\ j 
b + * V — ln U J - 2 ln(M£) 



X 



(V2«)'' 



-|7-Wt]-^(^ 



J! 



2pg 



-i-i 



£>7T 



sin|--W — 



2pg 



Equation ()135|) suggests that we define a new function ic- (£) by 



i*T n (x) = exp 



pln(p) +gln(g) /2po /g\ j 

i + e V T ln U J " 2 ln(w£) 



^ 12) (0- (136) 



Using (|13f)j) in (JIJ) we obtain, to leading order, 

(j + l) Rfl\ - v^,- + Rf-l = 
which has the independent solutions [HE] 

£f ) (0 = ~yDj (V^) , and ijf*> (0 = (±v^) (±i 

Using (|B3|) and matching with (|135|) we get 

pln(p) + q ln(g) 



iT n (x)~K ? (12) (0 = exp 



+ ^/^ln (J J - ^(pge) 



(t-H 2j t f-W¥ 



(135) 



exp(^) (137) 
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where the function A,- : 



is defined by 



'138) 



9.2.1 Matching the corner at (q, 1) and the interior of E 

Finally, we verify the matching between (|1U8|) and ()137j) . Using ()1 18|) and (|134|) in (|23|) we 

have, 

^ + g ln(g) + p ln(p) + piri 



111 I - I 7T1 

P 

2 



-1/2 



+| [1 - ki{pqje) - Tri] + v^j^i + y 



•0 g l n (o0 + P m (p) ~~ P 7 ™ 



In I-) Tri 

P 

2 



-1/2 



+| [1 - In(pgje) + Tri] - ^/2#i + y • 



Similarly, from (|2b|) we get 



1 



1 



and therefore 



z) ~ 



: exp 



gln(g) + p\n(p) 



+ ^V^ln(^^/ 2 + 



(139) 



x exp ^ — [1 — \&{pqje)\ \ cos 



2# 



J7T 

2 



Using the formula 

D n (0 ~ v^exp || pn(n) - 1]} cos [Vn£ - ^ 

in ([ISID yields 

pln(p) + gln(g) 



n — > oo 



^ 12) (0~exp 



^/^ln (jj - |ln(pge) 



, J 



^expj^[l-ln (j )] 
(140) 



cos 



pn 

e 



2pq 

e 



\Tj£- y) cos ( — -^\f— ) - sin f y^C - ? ) sil > i — - 



P7T 



2pg 

£ 



where we have used 



A, (0 



^ exp <j | [1 - ln(j)] j> sin ^ Vj£ } 



jit 



Equations f)139j) and ()140|) are identical, after regrouping terms. 
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Figure 2: A sketch of the different asymptotic regions. 

10 Summary and numerical results 

Below we summarize our results for the various asymptotic approximations to K n (x) as 



iV — > oo, with 



1 y z 

e = —, x=-, n = - 0<y,z<l. 

Nee 



[see Figure EI). 

1. Region I: n — 0(1), < y < 1, y 96 p. 



K n {x)~K£\y) = £ —{y-p) 



(see Figure El). 
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Figure 3: A comparison of K n (x) (solid curve) and Kn\y) (ooo) for n = 2 with e = 0.01 
and q = 0.64894783. 
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Figure 4: A comparison of K n (x) (solid curve) and K\ (rf) (ooo) for n = 2 with e = 0.01 
and g = 0.64894783. 



2. Region II: n = 0(1), y ~ p, y = p + r]y/2pqe, i] = 0(1). 

K n (x) ~ JtJOfo) = ^ (^) n/2 tf n fo) 



n! V 2 / 

where if n (rj) is the Hermite polynomial (see Figure HJ 
3. Region III: < y < K~(2:), < z < p, where 

Y ± (z)=p+(q-p)z±2zU , U (z) 



pq(l - z) 



K n (x) ~ K (y,z) = e 1/2 — = exp [e 1 ip (y, z)] L (y,z), 

v In 
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Figure 5: A comparison of K n (x) (solid curve) and K (y, z) (ooo) for n — 10 with e = 0.01 
and q = 0.34894783. 



with 
and 



l(U± ~p)(U ± + q) 
z [(U±) 2 - U$] 



(see Figure EJ). 



1 fp-y 



2\z 



+ q-p)±U(?-JL + q -p) l((o) 



4. Region IV: Y + (z) < y < 1, < z < q. 



K n {x) ~ K+(y, z) = e 1/2 ^= exp [e^iy, z)] L + {y 
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Figure 6: A comparison of K n (x) (solid curve) and K + (y, z) (ooo) for n — 10 with e = 0.01 
and q = 0.34894783. 



(see Figure IHI). 
5. Region V: x = 0(1), p < z < 1. 

K n (x) ~iT (5) (x,z) = 



,1/2 



271^^(1 - z) 



COS! TTX 



z — p 

p 



exp 



(f>o( z ) 



e x 



-T(x) sin(7rx) 



qe 



z — p 



exp 



[z — 1) ln(g) + niz 



IT Z — p 

where 

(fio(z) = [z — 1) ln(l — z) — z \si{z) + z ln(— p) 
and T(x) is the Gamma function (see Figure |7j). 
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Figure 7: A comparison of K n (x) (solid curve) and (x, z) (ooo) for n = 80 with e = 0.01 
and q = 0.74894783. 
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Figure 8: A comparison of K n (x) (solid curve) and K^ 6 '(x, u) (ooo) for n = 25 with e = 0.01 
and q = 0.74894783. Note that with this choice of parameters, u = 0.024265. 



6. Region VI: x = 0(1), z ~ p, z = p — u^/pqe, u = 0(1). 

£ l/2 



K n (x) ~ K {G \x,u) 



^2wpq 




DM 



x exp 



nip — gln(g) u^/pqTxi—u^/pqhi(q) u 2 

e 4 



where D x (u) is the parabolic cylinder function (see Figure |HJ). 
7. Region VII: < y < F~(2r), p < z < 1. 



K n (x) ~ K (7 \y,z) = exp 
(see Figure 01). 



/Try 



cos ^— j K + (y,z) + 2isin ^— ) K~(y,z) 
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Figure 10: A comparison of K n (x) (solid curve) and K^ 8 '(/3, z) (ooo) for n = 10 with e = 0.01 
and q = 0.34894783. 



Region VIII: y « Y~(z), < z < p, y = Y~(z) - /fe 2/3 , (3 = 0(1). 



K n (x)~K^(f3,z) = e 1 / 3 



exp 



U -q 



Ai [e 2 / 3 /?] 



0-1/3 



where 



%{z) = zni + (z-l) In (U ) + K~(z) In (U - q) + [l - K^)] In (C/ + p) 

6(Z) = ^(Uo+p) (U -q) 
and Ai (•) is the Airy function (see Figure EH)- 
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9. Region IX: y (z), p < z < 1. 

K n {x) ~K i9) {(3,z) = e 1/3 exp 
1 iT 1 / 3 



e 1 ipo{z) + In 



x 



[A+(/3, z)Ai (tf 2/3 /3) + iA"(/3, z )Bi ($ 2/3 f3)] 



where d(z) = —Q(z) 



A ± (/5, z) = exp 



2th [Y~ (z) - /fe 2 / 3 ] ] ± i 



and Ai (•) , Bi (•) are the Airy functions (see Figure ITT]) . 

10. Region X: Y~(z) < y < Y + {z), < z < 1. 

K n (x) ~ K( 10 )(y, *) = *) + z) 

(see Figure H2))- 

11. Region XI: n^N, n — N — 0<y<l, y^sq. 



K n {x)~K? 1 \v)=[*: )(-pT-' 



N-j ( 1 



Ny ( y 

cos (Nny) 1 



Ny 
N-jJ \q 



sec 



Figure H3J) . 



12. Region XII: n m N, y m q, y = q + £y/2pqe, £ = 0(1) 

pln(p) + gln(g) 



K n (x) ~ ir] 12) (0 = exp 

Id, (v^e) cos 



+ ^/^ln (jj _iln(pge) 



exp 



X 



where the function A,- : 



e V e 



is defined by 



1 a (c\ ■ ( prr c l 2 ^ 



(see Figure [Til). 
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Figure 11: A comparison of K n (x) (solid curve) and K^(/3, z) (ooo) for n = 40 with e = 0.02 
and q = 0.74894783. 
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Figure 12: A comparison of K n (x) (solid curve) and K^ 10 \y, z) (ooo) for n = 40 with e = 0.02 
and q = 0.74894783. 
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